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Abstract 

We study existence and uniqueness of distributional solutions w to the ordinary differ- 
ential equation -^^a{x) ■ — ^^^7-^^ + P{x)^^-^^i^ = g{x) with discontinuous coefficients and 

right-hand side. For example, if a and w are non-smooth the product a ■ w" has no obvious 
meaning. When interpreted on the most general level of the hierarchy of distributional prod- 
ucts discussed in [101 Chapter II], it turns out that existence of a solution w forces it to be at 
least continuously differentiable. Curiously, the choice of the distributional product concept 
is thus incompatible with the possibility of having a discontinuous displacement function as 
a solution. We also give conditions for unique solvability. 

Key words: ordinary differential equations with discontinuous coefficients, distributional so- 
lutions, multiplication of distributions. 
AMS 2000 Subject Classification: 46F10;34A36. 

1 Equation of the Euler-Bernoulli beam 

We consider an Euler Bernoulli rod under a distributed transversal force g and axial force P. The 
differential equation of equilibrium for the displacement w is given in [T] in the form 

Here, E is the modulus of elasticity, / is the moment of inertia, and I is the length of the rod. In 
our analysis we will allow for nonconstant, x-depended, even discontinuous coefRcients / and P. 
When there is a discontinuity in / at some point xo the rod can be considered to consist of two 
different, but connected, parts, i.e., EI{x) = EIi + H{x — xq){EI2 — EIi) where Ii 7^ I2 are the 
corresponding moments of inertia respectively and H denotes the Heaviside function. 

Equation ^ has been studied in (12j . where the authors discuss possible jump discontinuities 
at xo in the displacement 

A ■= wixQ+) - w{xo-) (2) 
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(as well as in the rotation 6 :— u''(a;o+) ~ w'{xq~)), where a sufhx + or — in the function 
argument denotes the limit from the right or left. A solution Ansatz of the form w — wi + 
H {x — xq) {w2 — wi) is then used, where wi and W2 solve the equation to the left and to the right 
of xo- In course of justifying this could be called a solution A was being forced to vanish in order 
to avoid ill-defined products involving a Dirac delta. 

Here, we investigate the corresponding mathematical issues left open: first, we analyze the 
possibility to give a meaning to the notion of 'distributional solution' in the context of the distri- 
butional product hierarchy described in [10', Chapter II] (see also the Appendix for a brief review); 
second, we show that indeed A necessarily has to vanish then, which is consistent with the calcula- 
tions in [12] : more precisely, iiw were to have a jump discontinuity then the model product [a-w"], 
which is the most general in the distributional product hierarchy, can not exist. Thus, in order 
to allow for solutions with jump discontinuities in the displacement one is forced to go beyond 
intrinsic distributional products and use, e.g. algebras of generalized functions (cf. [HHn]). For 
example, there has been active research on such issues for hyperbolic partial differential equations 
with discontinuous coefficients, where in certain cases non-existence of distributional solutions has 
been proved (cf . [1 [3 1 [H E [9] ) . 

For notational simplification and structural clarity we put A — EIi, B = EI2 (hence A B) 
and 

a{x) ^ A+{B- A)H{x - xq) = AH{xo - x) + BH{x - xq). (3) 

Then the governing differential equation with boundary conditions for the Euler-BernouUi rod 
with jump discontinuities in the bending read 



' ')+P{x)^A^^g{x), xe[0,l] (4) 



dx'^ I dx^ 

u;(0) = 0; w(/)=0; 

Mechanically, a global condition of equilibrium is expressed by equality of the bending moments 

Ehw"{xQ-) = El2w"{xo+). (6) 
We may use the substitution u = w" to lower the order of equation ([4]) and boundary conditions 

m 

^[a{x)-u(x))+P{x)u{x)=g{x) x G [0, 1] (7) 

with 

w(0)=u(l) = 0. (8) 

Remark 1.1. Note that the above substitution is equivalent to imposing the additional boundary 
problem -^w{x) = u(x) with w(Qi) = w{l) = 0, which is uniquely solvabile once u is determined 
by (I7][S1). In the sequel we will thus only consider u. 

In equation ([7]) the product of the distributions a and u arises. There are several concepts of 
partialy defined products in the space of distributions. In the current paper we use the so-called 
model product (cf. TO]) to give a meaning to the differential equation. 

Remark 1.2. (Comparison with L'^-operator theory.) The above boundary value problem 
([7]|S1) can as well be investigated in the classical functional analytic context of unbounded operators 
on i^([0, 1]). Singularities of the coefficient functions then have a significant influence on choices 
for an appropriate domain. In course of the current paper, we follow an intrinsic distribution 
theoretic view, which allows for a wider class of solutions, right-hand sides in the differential 
equation, as well as variations in the solution concept itself. 
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To illustrate the situation in an unbounded operator approach we briefly sketch the con- 
structions for the case where a is given by ([3]) and P is a real constant. It is natural to im- 
plement the boundary conditions ([5]) into the domain of the operator. Furthermore, we have 
to specify the meaning of the formal expression (au)" . Note that requiring that u belongs to 
the Sobolev space -ff^(]0, 1[) makes u" well-defined in L'^([0, 1]) and gives sense to the bound- 
ary conditions u(0) = = 0. Observe that under these hypotheses au — Au^ -I- Bu+, 
where m_ (resp. w+) vanishes to the right (resp. left) of xq and is continuously differentiable 
on the left (resp. right) up to xq. Thus, by Schwartz' formula ([TTl Chapitre II, §2]), we have 
{au)" — au" + [Bu'j^{xo) — Au'_{xq)) ■ S^,^ + {Bu+{xq) ~ Au-{xq)) ■ 6'^^, which is in I? only for u 
such that the coefficients of 5xq and 5'^^^ vanish. 

Therefore, we define the operator Tu := a ■ u" with domain 

D{T) := {u e i?^(]0, 1[) : u(0) = u{l) = 0,Bu^{xo) - Au+{xq) = 0,Bu'^{xq) - Au'_{xo) = 0}. 

It is straightforward to check that T is symmetric, i.e., D{T) C D{T*) and T* \d{t)= T, where T* 
denotes the adjoint of T. In fact, one can prove that T is self-adjoint along the following lines: Let 
(p G C^(]0, l[)r)D{T) and v e D{T*); interpreting i^-inner products (.|.) in terms of distributional 
actions (. , .) and vice versa we obtain {T*v, ip) = {(p\T*v) = (Tip\v) = {aip"\v) = ((aw)", (p), which 
implies that [av)" € L^, forcing that v belongs to and satisfies the conditions appearing 
in D{T) at Xq. Furthermore, integration by parts is then applicable with u € D(T) yielding 
{Tu\v) = {u\a ■ v") -\- Bu'{l)v{l) — Au'{0)v{0); since u i-^ {Tu\v) has to be a continuous linear 
functional with respect to the L^-norm v(0) and v{l) have to vanish. Hence v is in D{T) and 
T* = T. 

We observe that the original differential operator in Equation ([7]) is of the form T + PI , 
where / denotes the identity operator. Therefore, questions concerning uniqueness and existence 
of solutions to ([7][5]) when g G directly relate to spectral properties of T. One can view 
corresponding results obtained in Section 3 below in this context. 

2 Solution concept based on the model product 

We analyze the properties of a distributional solution u to problem ([7][8l) in detail when the product 
a • M is interpreted as a 'model product'. Throughout this and the following two sections we focus 
on regularity issues stemming from the highest order terms in the equation. Therefore we make 
the assumption that 

P is constant. 

We will remove this assumption and generalize our results in a final section allowing for jump 
discontinuities in P as well. 

Definition 2.1. Let I>'([0, 1]) {v e P'(R);suppw C [0, 1]}. We call u G i>'([0, 1]) a solution 
to (013^ */ the following holds: 

(Al) The model product [a ■ u] of u and a (defined as in JlOf, see also the Appendix) exists in 
(A2) The equation 

{[a ■ u])" + Pu^g (9) 

holds in V'{R). 

Remark 2.2. (i) The boundary conditions ([8]) are implemented into the definition of the space of 
prospective solutions T>'{[0, 1]) in the following sense: if u G T>'{[Q, 1]) happens to be a continuous 
function then u{0) = u{l) — 0. 

(ii) Note that (Al) is equivalent to the existence of the model products • u] and ■ u] where 
H-{x) — H{xq — x) and Hj^ — H{x — xq). 
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Lemma 2.3. (i) Let u G I>'([0, 1]) satisfy (Al-2) then [H^ ■ u] and [H+ ■ u] belong to i>'([0, 1]). 

(a) • (jip''] and • 5^xq] exist if and only if k — 0, in which case we have • Sxg] = 
[H^ ■ 6xa] = ■ (Cf. similar investigations in '5, Lemma 4]) 

Proof. Let ip e I?(R) with and tps{x) ip{x/e)/e be a model delta net (cf. [ID],(7.9)). 

(i) By definition [H± ■ u] = lim£^o(-ff± * Ve) ■ {u * ^e)- Let := {H± * (p^) ■ (u * (pe) and 
■0 £ I?(]R) with supp^/i n [0, 1] = then {w^, V") = {H± * (fe, * 'fie) ■ i')- Since supp(M * (p^) C 
[0, 1] + suppiy9£ C [— c?e, 1 + dg] for some d^ ^ (as £ ^ 0) we have {u * ip^) •■0 = and thus 

{We,(p) = 0. 

(ii) For any ip e I?(M) 

{[H^ ■ <5«], V) = lim((i/(.To - x) * ^e)(<5W * v'e), V') 



= lim 

£^0 e 



k+i 



<^(t)<^W(^^)^(a;)dtda; 



lim — r 

£^0 £* 



cannot be convergent for all as e ^ if fc ^ 0. In case /c = we obtain the formula [H- ■ 5xa\ — 
—Sxal^ by dominated convergence and the fact that J^J^ (p{t)ip{z) dzdt — —1/2. The proof for 
■ (5^0 ] is similar. ■ 

Theorem 2.4. Let u G 2?'([0, 1]) 6e a solution in the sense of Definition \2.i\ Then u is a locally 
integrable function. 



Proof. Step 1: Putting u_ 



" l(0,xa) 



and U4 



\{xo.l) 



yields 



Au_ +Pu^ = g |(o,a;o) 
Bu'l + Pu+ = g . 

Solving these two differential equations with constant coefficients we get 

U- — u-fi + u-p and u+ — u+h + u+p, 

where 



(10) 
(11) 



(12) 



and 



U-p{x) 



with a similar formula for uj^p(x) replacing P/A by P/B and integration limits from 1 to x. 
Here, u+h are smooth and u-p, il+p are absolutely continuous. Therefore U- and M4. are 

absolutely continuous functions on open subintervals (0,a;o) and (xq, 1) respectively. Also, by 
explicit formula, we see that u_(xo— ) :~ lima:^^;;,- u_(a;) and u+{xo+) :~ limj:^^;^^- M+(a;) exist. 
Step 2: Define u € L\^^{R) by 



u{x) — < 



-00 < X <Q 

u-{x) < X < xq 

U+{x) Xq < X < 1 

1 < x < 00 



(13) 
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We have that u G I''([0, 1]) and (m — u) \k\{xo}= 0. Therefore supp(u — u) = {xq}, which imphes 
that 

N 

u = u + Y,CkSi''J, CfceC,iVeNo. (14) 

fc=0 

By Lemma EH and Assumption (Al) iV = in Hence 

u^u + coS^g. (15) 

Step 3: By Assumption (A2) we now obtain 

{u ■ a)" = g-Pu^g~Pu- cP5^„ , (16) 

where g — Pu e Lj'^j,(R). Let w be a primitive function for g — Pu. Then w — cqPH(x — xq) is 
one for [u ■ a)' . Therefore 

u ■ a = W — cqP(x — xq)+, 

( X 2; > 

where xa- denote kink function, i.e. a;+ = < _ ^ „ and W is primitive function for w. Since 

W £ and the kink function is absolutely continuous we have that u • a is absolutely continuous. 
But then ^ and ^ imply that 

u • a = u • a + "^^xo {B — A) 

which is absolutely continuous if and only if cq = 0. This in turn yields u — u and therefore u is 
locally integrable as u is. ■ 

3 Existence and uniqueness of an L^([0, l])-solution 

As we have seen in the previous section, a distributional solution in the sense of Definition 12.11 
necessarily is a locally integrable function. In this case, we can interpret the product u • a as a 
duality product (cf. [TO] or the Appendix). We analyze this situation more closely. 

Proposition 3.1. If u G i^([0,l]) is a solution to (QQ) then u G C^([0, l]\{a;o}) and u has a 
jump at X — xq. 

Proof. If u e L^([0,1]) then the differential equation ([7]) yields (u • a)" e L^([0,1]), hence 
(u-a)' € Cabs([0, 1]) and thus u-a £ ^'^([0, 1])- Therefore we also have that u-a |[o,2:o)= l[o,xo)^ 
C^([0,a;o)) in turn u G C^{[0,xo)). Similarly, u € C"'^((a;o, 1]). Furthermore \imx-txo- ^ ' O'ix) — 
A ■ u{xq) and therefore u{xo—) = limj;^^;^- u(a:) exists. Similarly for u{xo+) — Vnnx^xo+ "fJ-ix) . 
But a • u is continuous, so that limj;^^;,,^ u ■ a{x) ~ lim2;_^2;o+ ' aix) and thus 

Au{xq-) = Bu(xqV), (17) 

which implies the global equilibrium condition If u G C([0, 1]) then p?| implies A = which 
contradicts the assumption I\ 7^/2- This means that u has to be discontinuous at xg. ■ 

Remark 3.2. As a matter of fact we have v! £ Cabs([0, l]\{a::o}). Indeed, since {u-a}' G Cabs([0, 1]) 
reasoning as above we obtain that u' is absolutely continuous off xo, so that u'{xo~), u'(xq+) exist 
and obtain 

Au'ixo-) = Bu'{xa+). (18) 

Now we are in a position to construct a solution to ([7][8l) . 

Lemma 3.3. For any choice of A> 0, B > 0, and < xq < 1 there exists a strictly increasing 
sequence {Pi)i^fi of positive real numbers Pi such that the following holds: 
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(i) IfP<0 
or 

(ii) if P > and P Pi for all I e 'N 

then there is a unique solution to ilO\) and II 11]) with w_(0) = and "+(1) — 0, which satisfies the 
stability conditions 



Au-{xo) = Bu+(xq) 
Au'_{xo) = Bu\{xo). 



(19) 
(20) 



Remark 3.4. In case P = Pi for some I E N the solution is not unique or even may fail to exist. 
Investigation of these cases seems possible in a direct way without requiring further analytical 
tools. 

Proof. Any solutions to pT7)l and pT|) are given by (H^). 

Case P < 0/ The solution formulae (fT2l) . adapted to the boundary conditions at and 1, give 



l-ix) ~ 2Ci sinh ^ —PjAx + u^^ix) 



and 
where 



u+(a;) = 27:>ieV--P/-B^ sinh ^J-P/B {x - 1) + u+p{x), 
1 



U-p{x) 



V^p/aJo 



;(t) sinh y/-P/A{x - r) dr 



and similarly to u+p{x) (replacing P/A with P/ B and integration limits from 1 to x). The stability 
conditions (|19II20I1 are equivalent to the linear system Hy — z with 



H 



2 A sinh ^^-P/Axo 



-2BeV--P/'B sinh ^-P/B{xo - 1) 



and 



2A^y-P/A cosh ^-P/Axa -2B^/^P/BeV^^ cosh ^-P/B{xo - 1 ) 



Bu+p{xo) - Au^p{xo) 
Bu'_^p{xo) ~ Au'_p{xo) 



(21) 









y := 




z := 



Further we have 



det H = 4:ABey^^ (^-^-P/B sinh ^-P/Ax^ cosh ^-P/B{xo - 1) 
+ y/-P/Asmh y/-P/B{xo - 1) cosh y/-P/Axo , 



Since 
and 



sinh ^-A/Pxo cosh ^-B/P^xq - 1) > 



sinh ^y-B/P{xo - 1) cosh ^^-A/Pxa < 



we have det H < hence unique solvability of the above linear system. 

Case P > 0: If P > then the solutions m_, u+ involve sin and cos (instead of sinh and cosh) 
and the determinant of the corresponding linear system Hy = z reads 



deti? = AABe^/^^ [-^JP/Bsm ^JP/Axq cos y/P/B{xo - 1) 
+ V'P7Isin ^/P/B{xo - 1) cos ^/P/Axa , 
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When det H we have the same situation as in the case P < 0. Observe that the set Zq of 
values for P such that any cosine factor occurring in the above determinant vanishes is at most 
countable. Apart from these values, to find P > for which det H — is equivalent to solving 

h{s) :— tan(s) + j///tan(/^s) = 0, 



where s = ^P/Axq > 0, — ^/AJB{1/xo — 1) > 0, and = xo/(l — xq) > 0. One observes 
that there is a countable discrete set of singularities of h, at which the limits from the left and 
right are +oo and — oo respectively. Since h is continuous otherwise, there is a countable set Zi of 
(positive) zeroes. To summarize, the union ZqU Zi makes up a sequence (P;)zgN with the required 
property. ■ 

Remark 3.5. We point out that the above proof of Lemma [3.31 does not give the minimum set 
of values P; to be removed. In fact, only those elements in Zq have to occur in (P;) which make 
both cosine factors vanish. Note that the latter can only happen, when ^y B / Axq / (1 — xq) is a 
rational number of the form {21 + l)/{2k + 1) with integers k, I. 

Theorem 3.6. Let P satisfy one of the conditions (i)-(ii) in Lemma \3.3l Let U- and m+ he the 
solutions to ilO\) and ill]) obtained in Lemma \3.3\ and define 

u.{x) = ( "-i^) " ^ , u+ix) = I ~ ° ^ ^ !^'^?! . (22) 

^ ^ \ a;e[a;o,l] ' I u+ix) xe[xa,l] ^ ' 

Then 

u{x) — u^{x) + u-^-{x) (23) 

is the unique solution to in the sense of Definition \2.I[ belongs to L^{[0, 1]) and satisfies the 

boundary conditions in the classical sense. 

Proof. Since a{x) — AH{xo — x) + BH{x — xq) we have 

{u ■ a){x) = Au^{x) + Bu^{x) 

which we will differentiate twice. Recall ([IH Chapitre II, §2]) that if a function / is in Cabs([0, l]\{xo}), 
such that lim.x^xo- fix) = /(xq— ) and lmix^xa+ f i^) = fixo+) exist, then the distributional 
derivative ^/ satisfies 

^fix) = fix) + {f{xo+) - fixo-)) ■ <5.o, (24) 

where /' denotes the (class of) function(s) in i^([0, 1]) equal to the pointwise derivative of / 
almost everywhere in [0, 1] \ {xq}- 
Therefore 

— (u • a) = Au'_ + Bu'j^ + (Pu(a;o+) — j4u(xo— )) • 5xq 

and 



^-^(u • a) = Au"_ + Bu'l + (Bu'ixf)+) - Au'ixf)—)) ■ S^q 



+ {Buixo+) - Auixo-)) ■ S',^^. (25) 



By construction we have that 



Au"ix) = { -e[0,xo] ^^.(^)^f ^ ^ xe[0,xo] 

^ I xe[xo,l] I i-Pu + g)ix) x£[xo,l] 



Thus dni) and HH) imply that 



-^iu-a) = -Pu + g. 



Note that u is continuous near the boundaries x ~ and a; = 1, thus the conditions ([S]) follow by 
construction. 

For uniqueness, we first observe that any solution u has to be in Lj^^ by Theorem 1 2. 41 Further- 
more, due to Proposition [3Tl] it also has to satisfy the stability conditions (|17m8p . Hence Lemma 
implies uniqueness. ■ 
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4 Generalization to discontinuous axial force 



We extended the analysis of the previous section to investigate solvability of the same type of 
differential equation 

(f 

[a{x) ■ u{x)] + Pu{x) = g{x) 
with boundary condition u(0) = ""(1) = 0, where the force P now is a jump function of the form 



P^Pi+H{x-xo){Pi^P2) 

with real numbers Pi, P2. 

As with constant P we obtain that any solution to the differential equation (in a sense similar 
to Definition 12. ip necessarily is a locally integrable function and continuously diffcrcntiable off xq 
with a jump dX x = xq. 

Remark 4.1. Note that condition (Al) in Definition 12.11 implies that the model product \P ■ u] 
exists. Therefore we will now require w to be a solution to the differential equation in the sense of 
Definition 12.11 with ([9]) replaced by 

([a-u])" + [P-u]^g. 

Theorem 4.2. (i) Let u G f([0, 1]) he a solution in the sense of Remark \4-.1\ Then u is a locally 
integrable function. 

(a) If u G L\^^ is a solution then u G C^([0, l]\{xo}) and u has a jump discontinuity at x = Xq- 
Furthermore, equations |j7| j and il8\) hold. 



Proof. Step 1: As in the proof of Theorem 12.41 we can set U- = u |(o.xo) ^-iid u+ = u \ {xo,i) ■ 
Then we have 

Au!_ + = g I (o,:ro) (26) 

Bw:;+P2W+ =5|(xo,i) (27) 
Solving these two differential equations with constant coefficients we get 

u_ = + u-p and u+ = u^h + u+p (28) 

where 



u-hix) = CieV^^- + C2e--F^\ u+^ix) = D^eV^^^ + D^e-^F^^ 

and 

u.,{x) = -^== { rg{r)eV^^(^-r)dr- P g{r)e-V^i^-^)dr), 
Z^ — Pl/A Jo Jo 

with a similar formula for u+p(x) replacing Pi /A by P2/B and integration limits from 1 to 
X. Again, are smooth and M_p,'u+p are absolutely continuous. Therefore S_ and 

are absolutely continuous on open subintervals (0,xo) and {xq,^)- Also, there exist u^(xo~) — 
lim2,^2,o_ M_(x) and u+{xo+) = lim^^^o^ u+{x). 

Step 2: Precisely as in Step 2 of the proof of Theorem 12.41 we obtain 



u = u + cSxg. (29) 

Step 3: Equation ^ and P = Pi + H{x - xo){Pi - P2) leads to Pu = Pi2_ + P2U+ + 
§Sxo{P2 — Pi)- Since PiM_ + P2U+ G Ll^^{M.) we have that w — J{g — Piw_ — P2W+) is absolutely 
continuous and its primitive function W is . The differential equation (u ■ a)" — g — Pu implies 

u-a = W^-^(P2-Pi)(a;-a;o)+, 
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which is absolutely continuous. The same arguments as in Theorem [23] yield that c = and hence 
u = u is locally integrable. This proves (i). 

For part (ii) we may reason as in the proof of Proposition 13.11 we get that second part of 
theorem is valid. ■ 

The construction of a solution rests on the following lemma which corresponds to Lemma 13.31 

Lemma 4.3. For any choice of A > 0, B > 0, and < xq < 1 there exist one- dimensional 
submanifolds A4 and Af of such that the following holds: 

(i) IfPi<0 andP2 < 

or 

(ii) if Pi > 0, Pa > and (Pi, P2) ^ M 
or 

(Hi) if Pi > 0, P2 < 0, and (Pi,P2) ^7V 

then there is a unique solution to \20fl and {2'T\j with S-(0) = and 2+(l) = 0, which satisfies the 
stability conditions 



Au^{xq) — Bu^{xq) 
Au'_{xo) = Bu'_^{xo). 



(30) 
(31) 

Remark 4.4. Similarly as in Remark 13.41 for the cases where (Pi,P2) belongs to or A/" the 
solution is not unique or may fail to exist, explicit investigation of which could be carried out 
along the lines of the following proof. 

Proof. As in the proof of Lemma anv solution to (pS)) and (|?7|) is given by 

In case Pi < and P2 < the solution formulae (j28p with boundary conditions at and 1, 

and stability conditions (j30p and (|3ip lead to a linear (2x2) system Hy ~ z with y and z as in 

Lemma 13.31 and 



det H = 4APeV--P2/B (-^-Pj/Ssinh ^-PijAx^ cosh yJ-P2/B{xa - 1) 



V-^iMsinh yJ-P2/B{xo - 1) cosh ^J-Pi/Ax^ < 



Therefore we have a unique solution. 

Case Pi > 0, P2 > 0: detH now reads 



AciH = iABeVP^ [-^/ P2/ B sin ^ Pi/ Axo cos ^/P2/B{xq - 1) 
+ /fVlsin ^P2/B{xo - 1) cos ^Pi/Axo) • 

Whenever this is nonzero we have a unique solution. To see 
where it vanishes let s — yJPijAxQ, t = \/I^/B{xo — 1), 
v = Xq / {1 — Xq) and analyze the function 

/(s, t) — j/t sin s cost + s sin i cos s. 

By direct inspection one deduces that grad/ is nonzero 
when / = which yields that the zero set Ai' = {{s,t) G 
R^;/(s,t) = 0} is a one-dimensional submanifold of M^. 
We set M = {(Pi, P2) G (s, t) G M'}. 




Figure: A^' is a union of infinitely 
many closed concentric curves (plot 
for the case = 6). 
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Case Pi > 0, P2 < 0: we have 



detH = AABe"/^ (^V-Pa/Ssinh ^/-Pl/Axn cos ^JP2/B{xo - 1) 
+ v/PiMsin ^P2/B{xo - 1) cosh /P^xo^ 



and as above one can show that the zero set is a one-dimensional submanifold A/" of R^. In the 
complement dct H =/= 0, solution exists and is unique. ■ 

Theorem 4.5. Let Pi and P2 satisfy one of the conditions (i)-(iii) in Lemma \4-.3\ Let u_{x) and 
Uj^{x) he the solutions to i26]) and {21^ obtained in Lemma \4.3\ and define 

r x&[0,xo] , x&[0,xo] , . 

X e [0^0, 1] "+^''^ ~\u+ix) x€ [xo, 1] ■ ^-^^^ 

Then 

u{x) — u-{x) + u+{x) (33) 
is a solution to (Q[3i with P ^ Pi + H{x - xo){Pi~ P2). 

Proof. Proceeding as in the proof of Theorem l3.6l we arrive at ((25)) . By construction we have 



Au" = j -^1;^ + ^ ^^P'^?! and Bu''=\ J 

y xe[a;o,lJ + [ -P2U + g x e [xo,l\ 

Employing ^ and §^ we get ^(u • a) = -(Pi + (P2 - Pi)H{x - xo))u + g. m 



5 Approximation by regular izat ion 

In this section we investigate the possibility to approximate the solution to ([7]), ([8]) using some 
regularization of the coefficient a{x). Throughout this section we will assume that P is constant 
and such that the solution to Q, ^ is unique. 

Let Qu := (au)" + Pu and consider the equation Qu — g E L^([0, 1]). Suppose that is a 
smooth regularization of the jump disconuity in a such that g ^^([0, 1]) with limg^o ai"''' — a*-"''' 
uniformly on compact subsets of [0, 1] \ {xq} for derivative orders j = 0, 1, 2. Let Q^u := [au^)" + 
Pu. 

Proposition 5.1. Let g e L^{[Q, 1]). 

(i) If u^ G f ([0, 1]) denotes the solution to 

QeU = g, u(0) = 0,m(1) = 0, 

then Ue belongs to the space ^C^([0, 1]) of continuously diifferentiable functions whose derivates 
are loccaly integrahle. 

(ii) Let u be the solution to Qu ~ g. Then u^ —f u uniformly on compact subsets of [0, 1] \ {2^0}- 

Proof. (i): Note that a^QeUe = a^g is equivalent to the Sturm Liouville problem Lu^ :— 
{pu'^y + qu^ — a^g with p :— a^ and q := a^{a'l + P). If Q^u^ — a^g then (a^u^)' = Oeg — a^{a'l + 
P) e i^([0, 1]) since a^g S i^([0, 1]) and a^ is C^. Therefore a^Wg is absolutely continuous. Since 
tte is bounded from below away from zero this implies that u'^ is absolutely continuous as well. 
Thus u, e ^C2([0,l]). 

(ii): Let Ve '■= u^ — u. Then we obtain 

Q,ve^{[[a,-a)-u])" -.^ f,, t,,(0) = t-,(l) = 0. (34) 
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By assumption and using that — a is off xq we have that — s- in L^q^{[0, 1] \ {xq}). 
Integrating twice in gives 

{a,v,){x)= {x^y)f,{y)dy + P {x - y)v,dy + a,{0)MO) + {aeVe)' (0) x. (35) 
Jq Jo 

Let fiT be a compact subset of [0,xo). Since u and Ug are AC^ on if (as noted in Remark 13.21 
and the prooof above) and both beiong to 2?'([0, 1]) we have that Ve{0) — v'^{0) — 0. Furthermore, 
by \{aeV^){x)\ > A\v^\ > we obtain 

1+p r 1+p r 

\ve{x)\<^—J \x - y\\fe{y)\dy -\ —J \x - y\\v,{y)\dy. 

Applying Gronwall's inequahty we get 

\vAx)\ < ^ £ \x - y\\My)\dy < WfeWmK)- 

As noted above, we have ^ in Ll^^{[0, 1] \ {xq}), thus ^ uniformly on K. 

The reasoning in case of a compact subset contained in (xq, 1] is similar. Since an arbitrary 
compact subset of [0, 1] \ {xq} is the disjoint union of two compact subsets in either part of 
[0, 1] \ {xq} the assertion is proved. ■ 

We illustrate the convergence in an example: we put A ^ 1, B ~ 2, P ^ 1, xq = 1/2, and 
design as the function, which is defined by a fifth order odd polynomial in [— e, e] and 
equal to a otherwise. As right-hand side we choose g{x) = — cos(llx)/-y/|a; — 2/3|. Note that 
g G L^([0, 1]) \ i^([0, 1]). The following plots show the regularized coefficients and corresponding 
solutions for parameter values e = 1/10, 1/30, 1/100. 




A Appendix: Hierarchy of distributional products 

For convenience of the reader we briefly review the basic definitions of the coherent distributional 
products described in [TUl Chapter II] in terms of a hierarchy. All these products yield the classical 
multiplication when restricted to smooth functions. 

We use fl to denote an open subset of and u for the Fourier transform of u. 

The most elementary product in this context is C°° • V , the product of a smooth function and 
a distribution, defined as the adjoint of multiplication by a smooth function in the test function 
space. 

Disjoint singular support: Assume that u, v are in I?'(ri) with disjoint singular supports. 
Then for any x G Q there is a neighborhood Qx and a function /a; € ^{Qx) such that either f^u 
or fxV is smooth. Then in ^l^ the product of u and v can be defined in the sense of C°° • V and by 
the localization properties of V (cf. [3], subsect. 2.2) this consistently defines a distribution in H.. 
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Wave front set condition: Let u G V'{n) and (xo,Co) € T*n \ := {(x,^) | x G r2,^ 7^ 0} 
(the cotangent bundle over fl with the zero section removed), u is said to be microlocally regular 
at (xqi^o) if there is 1^9 G (/3(a;o) 7^ 0, and an open cone F with axial vector ^0 such that 

is rapidly decreasing in V. WF(w) is the closed subset of T*il \ where u is not microlocally 
regular. 

If u, f G I?'(r2) their wave front sets are said to be in favorable position if (x,^) G WF(u) 
implies that (x, — ^) ^ WF(u). In this case the product of u and v can be defined as the pullback 
of the tensor product u ® v ^ V {VI x. Vl) hy the diagonal map Vl ilxfi, x ^ {x.,x) (cf. [3], 
Thm. 8.2.10). 



Fourier product: Given two distributions u, v G 'D'{U) we say that their Fourier product exists 
if for every a; G there is an open neighborhood and fx G V^fl), fx = 1 on fix, such that the 
j5^'-convolution of f^u and f^v exists. Locally near x, the product of u and v is then defined to be 
the inverse Fourier transform of f^u * f^v (for a definition of =y'-convolvability see [lOj . sect. 6). 



Duality products: Let X be a normal space of distributions, that is I? C X C I?' and V is 
dense in X. Assume that the dual space X' is (equipped with a locally convex topology so that it 
becomes) normal as well and that multiplication with a fixed element in V induces a continuous 
linear map both from X into X and from X' into X'. 

For any normal space of distributions Y denote by lioc the set of distributions v E T>' such 
that ipv G Y for all t/' G 2^- If w G (X')ioc and v G X\oc then the product of u and v can be defined 
by 

{u ■ v,ip) := {xu,ipv) 

for ij) £ T) and x chosen arbitrarily with x = 1 on supp('0). Note that in the above definition 
the left hand side denotes a (V,!)) pairing while the right hand side uses the pairing {X',X). 



Strict and model products: The basic idea is to regularize one or both factors by convolution, 
perform the multiplication in the sense C°°-I?' or C°°-C°°, and try to take the limit. The regularizing 
convolutions are carried out with two principal types of moUifiers. 
A net {p')e>o in I?(K") is called strict delta net if 

supp(p'') ^ {0} as £ -> (36) 

J p^x) dx = 1 for aU £ > (37) 

J \p'^{x)\dx is bounded independently of £. (38) 
A model delta net is given by specifying ip G I?(R") with J ip{x) dx = 1 and defining {ip^)^yQ 

by tPeix) = V3(x/£)/£". 

Consider the following four possibilities to define a product of u and v: 

u ■ [v] = \imu{v*p^) (1) 

£—♦0 

[u] ■ V = lim(u * p'^)v (2) 

£—♦0 

[u] . [w] = \im(u * p^){v * a'^) (3) 

[u ■ v] — \im{u * p'^){v * p'^) (4) 

£—♦0 

where the limit is required to exist in I?'(i?") and independent of the choice of (p^)e>o and ((T')e>o 
in the class of strict, resp. model, delta nets. This defines four types of so-called strict, resp. model, 
products. Since the definitions (l)-(3) turn out to be equivalent when using strict, resp. model, 
delta nets (cf. fTO], Thms. 7.2 and 7.11) we distinguish only the following four products: strict 
product (l)-(3), strict product (4), model product (l)-(3), and model product (4). 
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Coherence properties: The 

various products satisfy coherence 
properties and can be brought 
into the foUowing hierarchy table. 
Here, an arrow indicates that a 
product definition is contained 
and consistent with its successor 
in the graph. All products shown 
generalize the multiplication 

c°° • v. 
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strict product (l)-(3) ► strict product (4) 
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model product (l)-(3) >- model product (4) 
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